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Abstract 

The usual method of solving the free particle Dirac equation results in the so called continuum 
energy solutions, with \E\ > mc^. Here, we take a different approach and find a set of solutions with 
quantized energies given by E = ±2(j + 1/2) mc^. 



Part 1: Solving the Dirac Equation as a Differential Equation 
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The free-particle Dirac equation 
Can be written in spherical coordinates as 

(^i h do — •y^ m c) ip 

where 



d .a' Li 

-in — + * 

or r 



(1) 



a 



4x4 



a. 



2x2 



0", 



2x2 



and (T2x2 3''^ the usual 2x2 Pauli matrices. 
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We demonstrate the transformation leading to ([T]) explicitly in Appendix A. By using the properties of 
the angular momentum operators we can write 



L3 L_ 

L+ -L3 

L3 L_ 

L+ -L3 
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(When the 2x2 submatrices on the diagonal are identical, and the off-diagonal elements are null, we will 
sometimes alternate between the 2x2 and 4x4 version, depending on context). The 4 x 4 e object is 
given by 

' A 
A 
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with the 2 X 2 A submatrices given by 



A 



cos 9 sinOe '^'^ 
sinOe^'^ —cos 9 



We will now construct an ansatz to solve ([T]), seeking solutions which are eigentates of H, J^, J^. 
First, consider the following 2-sp/V7or objects which can be constructed from spherical harmonics. The 
first is called 6'+: 

a/J - "^j yj-l/2,mj+l/2 _ 



(3) 



And a second object that we call 9 



0' + 1 - ?7i~ y j+l/2,m,-l/2 
-a/J + l+mj Fj+l/2,m,+l/2 



(4) 



We will see that these objects have total angular momentum j = 1/2,3/2.... The index nij is half- 
integral and runs from rrij = j, 1/2, —1/2, ... — j. (Notice that the indices of the spherical harmonics 
are integral as required). (For a reference, see Bjorken and Drell, pg. 53). 

These 2-spinors obey the relations 



a'L,9j = -h{j + 3/2)9] 



and 



a'L,9+= h{j -1/2)9+ 

Notice how the 6'+'^ behave slightly differently when acted upon by cr*Lj, although they are each separately 
eigenvectors of it - this can be verified explicity by using the representation of a*Lj shown in Eq. ([2]), 
which we demonstrate in Appendix B. 
They both have the same total J^: 

J^9+'- = {L^ + 3/^74 + ha' Li) 9+'- = h^j{j + 1) 9+'- 

They also both have the same value of = hrrij. Additionally they obey 
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That is, they are conjugate with respect to the A operator. These last two identities follow from the 
spherical harmonic recursion relations and are demonstrated in Appendix B. 

We find that we can construct two 4-sp/nor solution classes to the Dirac Equation ([T]) from these: 
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and 
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where the F^' (r) and (r) are c-ntvmber coefficients which are functions of radius only. 
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Let's investigate the first solution tpj and substitute into the Dirac equation ([T]). We get: 
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which becomes 



{Ee/c-mc)Ff{r)ej 
{Ejc + mc)GUr) 9' 



-i hd/dr - i h{j + 3/2) /r) Fr(r) Oj 
-i hd/dr + i h{j - 1 /2) /r) G+(r) 6+ 



{-I hd/dr + % h{j - l/2)/r) G+(r) ^J-^9J 
{-t hd/dr - i h{j + 3/2) /r) Fr [r) J^^Of 



We can rearrange the terms and separate the matrix equation into two relations between the coefficients 



of the el 



{d/dr - (j - 1/2) /r) G+ 
{d/dr + (j + 3/2) /r) F~ 



(Ee/c—mc) ■ 



-i h 
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Now, the spherical Bessel (or Hankel) functions obey the coupled recursion relations: 

{d/dr - {n - l)/r) jn-i = -njn 
{d/dr + {n + l)/r) j„ = nj^-i 

If we make the index substitution n = j + 1/2, then these recursion relations become 



(a/ar-(j-l/2)/r)j,_v2 = 
(a/ar + (j + 3/2)/r)j„+i/2 = 

So that we can solve the system given by Eq. ([8]) by identifying 



e Ji-i/2l 



Kr 
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(allowing for an arbitrary phase between F and G). This means that the wavenumber k must obey two 
relations: 

{EJc-mc) 
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which can be equated to each other to eliminate k, giving the following relation satisfied by the energy E, 

2J + l+jV25 
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For the solution tpj to be a steady state solution, the energy E^ must be real. This restricts 5 to either 
5 = n or 5 = n/2. In the first case 5 = vr we get: 



E. 



2 (j + 1/2) mc^ 



with corresponding k from ( IT2l) or ([T3 



2m c 

~h~ 



(15) 



(16) 



In the case 6 = tt/2 we get solutions with an inverse relationship between energy and angular momentum, 
which for now we will dismiss as unphysical. 

So by solving the Dirac equation in spherical coordinates, and requiring the solutions to be eigenstates 
of H, .P and J^, we find that the energy is quantized. By using the solutions 'ip'j, we find a similar 
quantization of energies, except now with negative energy values. We can summarize the solutions for 
both 'ipj and in one expression: 

E = ±2(j + l/2)mc^ (17) 

where the index j = 1/2, 3/2, 5/2... with a degeneracy of 2 j + 1 for each index j. 

Now we must address the discrepency between our solutions, and the traditional approach which results 
in the continuum of energy solutions. The beginning of the discrepancy comes in the choice of the definition 
of the 9^'^ 2-spinors, as in Eqs. [31 IH The traditional approach, as taken in Bjorken & Drell, [2J, Messiah, 

or just about any other textbook on advanced quantum mechanics is to define 9^'^ like this: 
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1 



- rrij +1 y j+l/2,m,-l/2 



Q+ ^ _l_ 



VJ + yj~l/2,m,^l/2 
a/J - '^j ^i-l/2,m,+l/2 



These additional prefactors of and result in 2-spinor objects which are seperately prenor- 

malized with respect to the angular variables: 

^ 9^9dn = 1 



In addition, instead of Eq.O the objects now obey a different conjugate relationship with respect to the 
A operator: 

\9+ = 9-, \9- = 9+ 

This in effect wasties away the the j dependant factors from the problem. Furthermore, instead of choosing 
a unitary phase factor e**^ between F and G, as in Eq. [TTl the traditional approach is to allow for a factor 
of arbitrary magnitude, call it v. So instead of Eqs. [T2]and[T31 we have: 

^ = l(E/c-mc) 
V in 

. = v^^l^±^ (19) 
—i n 

The j dependance is gone, and the factor v, related to the relative normalization between the upper two 
and lower two components introduces a third variable into a set of two equations with three unkowns. The 
procedure continues by choosing k to be an arbitrary parameter which can take on any continuous value 
from to ± oo , and then solving for E and v in terms of k. E is determined from: 

E{k) = ±v/(MVTm2^ 

A continuous function of k which has a minimum magnitude of mc^. (This relation between E and k is 
also satisfied by the discrete solutions, except that k is fixed at discrete values). 

So then, which approach is correct? One might proclaim that the traditional approach is to be favored 
because it is more in line with experimental fact. We know that free electrons do have continuous energies, 
as even the simplest cathode ray experiments seem to demonstrate. Furthermore, the existence of a 
spectrum of lieavy electrons with masses of 2m, 4m, etc. has presumably not been observed. But as we 
will see later, these heavy electron states will still have the same charge-to-mass ratio of e/m, so they would 
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be indistinguishable from regular electrons in experiments using electric and magnetic fields. Furthermore, 
note that this solution approach which results in discrete energies has been undertaken in the rest frame 
of the electron. The inclusion of Lorentz boosts will still allow for a continuum of possible momenta and 
total energy. 

To resolve the dillema, we must determine if the j-dependant factors which result from the conjugate 
relations between the 9^'^ as in Eq. |5] are inherently necessary and natural, or if a normalization which 
eliminates them is correct. To end the ambiguity, we will turn to a method of solving the Dirac equation 
other than solving it as a differential equation. 



Part 2: Solving the Free-Particle Dirac Equation using Heisenberg 
Operator Methods 

To begin, we note that the Dirac 7 matrices can be viewed as representations of the four unit vectors of 
space-time. That is to say: 

xo = 7o, xi = 7i, X2 = 72, ^^3 = 73 (20) 

These unit vectors obey: 

Xf,- x^ = -{-ff,-f^ + -f^-ff,) = g^^ (21) 
Thus a four-vector object itself amounts to a 4 x 4 matrix: 



A = A^x^ = A^'y^ 



(22) 



With this identification, the dot product and cross products between two vectors are represented by 

1 . r r ^ i 



A-B = ^{A,B}, Axi? = -7V 



A,B 



(23) 



The symmetric dot product of two vectors is proportional to their anticommutator, and the antisymmtric 
cross product of two vectors is related to their commutator. (An advantage of this representation of the 
unit vectors as matrices instead of as column vectors is that it allows for an algebraic formulation of the 
cross product which can be written on a two-dimensional piece of papei0). 

We see that the Dirac operator ifi'y'^d^ is the equivalent of the vector gradient operator ifiV . Thus 
the Dirac equation is the statement that the gradient of ^ is proportional to i}): 

ihV ip = mclp 

It's interesting to note that in this matrix representation of vectors, the existence of spinor columns and rows 
which the vector operators can act upon is implicit. In other words, quantum mechanics is an implication 
of the matrix algebra of vectors. 

Now, consider the four-vector which represents position: 

X = x^-i^ 



■"^For example, no algebraic operations beginning with 
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can result in 
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The four-velocity can be calculated from the first derivative of this with respect to proper time: 

d 



X = —X = r + C7 
dr 

(In the rest frame of the electron, ^ = ^). Now, note that 
derivative with respect to time: 

= — T'y^c — T'y^€ 

We must now use the Heisenberg equation of motion to calculate r and e 



r'j'^e. So again, taking the 

(24) 

(25) 



From Eq. ([T]) the Hamiltonian can be written: 



H = J mc + c e —ih— — \- i 
\ or 

Remembering that r is a diagonal scalar, we can calculate r as 

+ ce (^-th^ + i^^ ,r 



ih 



n 2 
7 mc 



(26) 



-ce 



Now, to calculate e, we keep in mind that e is a full matrix operator, but a function of the angular variables 
only: 













r ' 



ih 



7° mc^ + c e 



(27) 



We can now substitute the results of (!26|) and (1271) into (!24|) . Using 7°e = — €7°, = 1 and 7°e — €7" 
27°e, we obtain after some manipulation: 

I r ?■ r 1 '2,1710 

r = -cY + -7*^ [e, e a'Li\ —re 

Remembering x = f+ 07° we get an expression for the four-velocity vector operator.: 



X = [e, ea Li\ —re 



ih 



(28) 



We will be more interested in the time-component of this four vector, which we can obtain from: 
After making use of facts like {e, 7°} = and [7°,cr*Lj] = 0, we obtain after some manipulation: 



(29) 



Now let's investigate the action of the xq operator on the spinor solutions of the Dirac equation. First, 
consider ijjj, (suppressing time dependence of ijj): 



[a' Li - ta'Li t) 



F-e- 



(30) 



A 
A 

in the solutions of Part 1. That is: 



Consider the operator e 



It is the action of A on the 9 + which is the source of the ambiguity 



s 



Ar = s9+,X9^ 

In the traditional approach, the factor s = 1. In Part 1, we chose to let s 



(31) 

For now, let us leave 



s unspecified, and carry thru the calculation in ([30 



-h{j + 3/2)F-9- 
h{j - l/2)G+9+ 
-h{j + 3/2)F-9- " 
h{] - l/2)G+9+ 
-h{j + 3/2)F-9 
h{] - l/2)G+r 
-2Ti{] + 1/2) " 

2h{] + 1/2) 



l/sG+9 
sF'9+ 
{l/s)n{j + 3/2)G+9 
sh{j - l/2)F-9+ 
h{j-l/2)F-9- 
-h{j + 3/2)G+9 
F- 9- 
G+ 9+ 



(32) 



-2(j + l/2)c7o^| 

Notice how the result is independent of the factor s. The double application of the operator e in (I3U1) 
results in factors of ^ x s causing the s dependence to disappear. However, since the j factors appear 
in a manner independant of the value of s, we can say that the presence of the j factors is natural and 
inherent. 



Now let's calculate the action of xq on ip^. 



Xq 



^{a'Li-ea'Li e) 



h{j - l/2)F+9^ 
-h{j + 3/2)G-9 
h{j-l/2)F+9+ 
-h{j + 3/2)G-9- 
h{j - l/2)F+9^ 
-h{] + 3/2)G-^ 
2?i{3 + 1/2) 




F+9^ 
G-9- 







sG-9+ 
{l/s)F+9 
shQ - l/2)G-9+' 
-{l/s)hU + 3/2)F+9 
-h{j + 3/2)F+9 

h{j - i/2)G-r 

F+9+ 
G- 9- 



(33) 



-2nij + 1/2) 

2(j + l/2)c7o^J 

Again, the s factors cancel out, and now we have a result which is proportional to j + 1/2. In general we 
have: 



i:o^;'^ = T2(j + l/2)c7oV; 



(34) 
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That is, for the electron states ^| the negative sign pertains, while for the positron states, the sign is 
positive. We would almost be ready to interpret Xq as the operator E/mc if it were not for the fact that 
the overall sign is the opposite of that in Eq. (IlTl) . We will see in a moment that this should be so, and 
is due to the fact that electrons have positive energy, but negative charge, and positrons have negative 
energy but positive charge. Apparently, the correct identification of the E operator, to account for the 
minus sign is 



E = -mc^^xo = — — 7° (e a^Li e - (T^Li) 
h ^ ' 

This operator is the Heisenberg counterpart to the operator H = ihd/dt. 

Let's see how this plays out further. Consider that the classical current is defined as 

Jfj, — P ^fj, 

We can then see that there is a correspondence with this classical definition if we premultiply Eq. (!34f) with 
the electron charge e and tp: 



e,p 



pxo = er'^'ior'" = Te2{j + 1/2) ^'^loi^ 

If we have already found a suitable normalization of ^, (which will discuss in Part 3), then the term on the 
far right, which we know to be the normalized Dirac current jo = ip'joip, underscores this correspondence, 
and allows us to discern that the charges of these states are quantized as g = =Fe2(j + 1/2). 
Now, let's again premultiply Eq. ( 1341) . but this time with 'ip'^ instead of ^p^. 

r^r^^x^rf = T2(j + l/2)c^^'Pt7oV'l'^ (35) 

The term on the far right is the invariant ^ip, so therefore, the term on the left, which is the expectation 
value of Xq in the rest frame of the electron, must also be an invariant. 

<xo>= I r'"^ xo = T 2(j + l/2)c [ r^"^ 70 (36) 



(For now, we have left the metric and bounds of the integration unspecified - we will deal with this in Part 
3). If we can find a way to normalize the solutions ip^'^ such that: 

Then we will have shown that 

< io >e,p= T 2(j + l/2)c 

Apparently, the correct correspondence principle between the energy and the expectation value of the zeroth 
component of four-velocity is 

(The presence of the minus sign would not be noticed in classical calculations, since they always involve 
the square of the energy or zeroth component of four-velocity). Thus the j dependent terms which we 
discovered in Part 1 are natural, and prenormalizing the 9^'^ spinors to eliminate j is inconsistent with 
this Heisenberg operator approach. 
Thus the relation 

Ee,p = ±2{j + l/2)moc' 
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stands, and the ambiguity mentioned in Part 1 is resolved. We can interpret this as the prediction of a 
spectrum of heavy electrons and positrons, with inertia! masses of m = 2{j + 1/2) mo- 

Note that since the charge of the state, q, is quantized in the same way as the rest-mass, the charge- 
to-mass ratio is constant: 

q _ T2(j + l/2)e _^ 

m ~ 2(j + l/2)mo "^^^""^ 

Which is consistent with the experimental fact that all electrons or positrons have the same charge-to- mass 
ratio. 



Solution Table 



= 2{j + l/2)mc^, Ep = -2{j + l/2)mc^ 



2mc 



jj+i/2{'^jr)0j 
i3j-i/2{Kjr)e^ 



jj-i/2{njr)9J 



- VJ +mj + l rj+l/2,m,+l/2 



+ rUj lj-l/2,mj-l/2 
. Vj ~ »^j^j-l/2,m3+l/2 



Table 1: A table summarizing all the solutions of the free particle Dirac equation. For each j value, there 
is a degeneracy of 2 ']+! . 



Part 3: The Normalization of the Solutions over a Discrete Lattice 

Now the volume integration of the spherical Bessel functions is divergent: 

jf{K,r)r^dr — oo 



/■ 



So that the solutions to the free particle Dirac equation are not normalizable in the usual way. Both the 
traditional continuum solutions, and the solutions we presented in Part 1 have this feature. We are aware of 
the usual techniques used to deal with this divergence, such as limiting the extent of the volume integration, 
and then taking limits of ratios, etc., but our goal here is to introduce a new method of normalizing the 
free particle solutions. 

All of the spherical Bessel functions have well defined, countable sets of zeroes. For example, jo{Kr) = 
sii^, is zero whenever r = —, n = 1,2,3.... Higher order Bessel functions like ji(Kr) = - ^2M£ 
have sets of zeroes for which there aren't simple formulae like for jq, but they are nonetheless also discrete 
and countably infinite. 

Interestingly, although no simple formula for the zeros of ji exists, there is a simple result when ji is 
evaluated at the zeroes of jo- That is: 

, nvr, , , sinnvr cosnTr , 1 

ji 0^ — = = = -ir — 



This allows us to conclude: 



n=l n=l 
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In addition it is also true that 



n=l 



since each term in the sum is zero. 

Now consider the j — 1/2, positive energy, negative charge, electron solution: 



j=l/2 



Noting that /^^ 9j^9j dQ = 2j + 2 and that J^^ ofo^ dQ = 2j, we see that 
Similarly, for the the j — 1/2 , negative energy, positive charge, positron solution: 



d/2 = 



We get that: 



So, both results give 1/2: 



E / ^?/2%rf^^ = 0.(l) + i-(l + 2) = i 



We can also calculate the value ip ■ ip when summed over the lattice of nn. For the electron state '4'i/2' 
we get 

E / ri/2ri/2d^=Yl I ^V7Vi^/2c?^^ = ^-(i+2)-o-(i) = i 

nTT t/ 47r 

But for the positron state V'1/2 we get: 

^^nJL ^47r ^ 

The opposite sign, but same magnitude! Apparently, the invariant quantity ^ - ^ is related to the (negative 
of the) invariant charge of the state. 

Now, so far, the fact that these sums converge for j = 1/2 is interesting, but this approach is not really 
useful as a normalization scheme, unless we can prove that all solutions -0^ can be normalized in this way. 
May we find a formula: 



E 

zeroes of jo 



^fp^'/dn = dj 



4tt 



(37) 



That is, can we prove that the sums converge to a finite coefficient for all the solutions of the free particle 
Dirac equation? Discovering a general formula for the dj would be nice, but not necessary. We will see 
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that the recursion formula for the Bessel functions will allow us to discover any dj by building up from 
di/2. Expanding the above formula a little more explicitly we have for electrons 



oo 



d^^ = (2j + 2) ■ 5^ j|+v2(^^) - (2J) ■ (38) 



n=l n=l 



and for positrons: 



oo 



< = (2J) ■ - (2J + 2) ■ (39) 

n=l n=l 

which makes it obvious that rf^ = —dj. We also see that to prove that the dj converge we must investigate 
the properties of E^r=ii|±i/2(^^)- 

Consider that the spherical Bessel functions obey the following recursion formula: 

jn{x) = ^ ^ jn-l(x) - 3n-2{x) (40) 

This allows us to relate the value of any j„ back to the values of ji and jo. by iteratively applying the 
above recursion relation. If we choose to evaluate jVi at one of the zeroes of jo. the result becomes even 
simpler. For example, for j2in7i) we have: 

■ f ^ 2-2-1 .... ,^0 3 (-1)"+^ (-l)"+^3 
Cleary, j2 can also be normalized over the integer lattice defined by the zeores of jo: 

32 i^TT) _ 2^ _ ■ _ 

n=l n=l 

Just glancing at Eq. |40] allows us to infer that since Xl^iiK^^) ^""^ l^^th converge, any 

sum involving higher j„ must also converge, since the power of mx in the denominator is always higher. 
But can we discover a general formula? Perhaps, and I have not yet done so, but I don't need to to make 
an intriguing observation. We have calculated (ii/2 = 1/2. Let's continue on and calculate ^3/2: 

00 o °° -1-1 

dl,2 = (2 ■ I + 2) . ^ jKnvr) - (2 ■ |) ■ j\^(n7r) = 5 • - 3 ■ - = 

n=l n=l 

An interesting coincidence this is, how the terms involving the Bessel sums interacted with the terms 
resulting from the integration of the angular variables to result in exactly zero. If we continue on in this 
manner, we find that 6/5/2, ^7/2. '^9/2 are all also zero. The question that arises then is are all the dj other 
than di/2 zero?. Since I don't have a general formula for the Bessel sums, I have to turn to numerical 
methods to answer this question. And the answer I have found is: almost all of the dj are zero, except 
for a very sparse subset. The next non-zero dj after di/2 occurs when j = 498|, followed immediately by 
j = 499|. All the other dj in between, at least as far as the Mathematica routines tell me, are exactly 
zero! Apparently, these states lying in between j = 1/2 and j = 498| have an invariant norm which, 
defined as in Eq. |371 is zero. An implication is that, since these states have zero norm, we might expect 
their existence to be suppressed. Continuing beyond j = 499^, the dj are again entirely zero, until we get 
to j = 574i and then j = 575^. See Table [2] for a list of the first few non-zero dj. 
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3 




1/2 


1/2 


498 1/2 


166.38 


/ion 1 /o 
4yy 1/ z 


— iDD.OO 


5741/2 


1.05071 X 10^36 


5751/2 


-1.05071 X 


6531/2 


1.59431 X 10^4° 


6541/2 


-1.59431 X 10^40 



Table 2: A table of the first few non-zero dj. 



If there is any validity to this normalization scheme, then these results seem to imply that most of the 
particle solutions E = ±2(j + l/2)mc^ are suppressed, except for a very sparse spectrum with non-zero 
norms dj. Now, we know from the results of particle physics that a sparse spectrum of particles does 
indeed exist, but I have not yet been able to match up the known particles such as the muons, pions, etc. 
with these dj. However, I have still presented this method hoping that it, or a slight modification of it, 
may be able to predict the known experimental spectrum of particles. I know that this is not the complete 
answer yet, because it results only from the free particle Dirac equation in the absence of any potentials. 
In a future work I hope to consider the effect of different types of potentials on the solution spectrum 
and the resulting dj. 



Conclusion 

Now, the result of the bound-state Cou/omb-potential Dirac equation is that the internal energies of the 
Hydrogen atom are quantized as En ~ mc^ — ^^'^i^ ■ Thus, the energy of bound-state electrons are very 

nearly E = mc^, for even the most deeply bound n = 1 electrons. I have found that the factors of y'' ^y- 
which resulted in the ambiquity mentioned in Part 1, do not come into play in the Hydrogen atom, and so 
my conclusion is that the traditional solutions to the Hydrogen atom problem are correct, as they must be 
(excluding hyperfine considerations). In a future work, I will apply the Heisenberg operator method to the 
problem of the Hydrogen atom, and show that it gives the same result as directly solving the differential 
equation, as I did in Part 2 of this work. We will see that since the covariant operator H undergoes 
a transformation H — > ih-^ — ^ as a result of the Coulomb potential, there must be a corresponding 
transformation of the contravariant Heisenberg operator x^. Discovering the laws of transformation of the 
Heisenberg operators x which result from the transformation of the covariant ihd^ + will be 

one central goal of a subsequent article. 

Finally, the results from Parts 1 and 2 show that the lowest energy of a free particle state is = 
2(1/2 + l/2)mc^ = 2mc^ , so that the first excited state free electron, has an energy of mc^ greater than 
a typical bound electron: 2mc^ — mc^ = mc^. But experimental work with cathode rays, the photoelectric 
effect and other elementary experiments seem to imply that free electrons have an energy always greater 
than mc^, say Ef^ee — mc^ + AE, where AE can range from a few eV, to several hundred-thousands of eV, 
at which point the free electrons begin to emit 7-rays and cause pair production. We might instead consider 
these free electrons to have energies of £^ = 2mc^ — AE' . Where AE' results from the perturbation of the 
free-particle state by a potential field. Again, we will address this in more detail in a subsequent article. 
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Appendix A 



To demonstrate the transformations leading to Eq. ([T]), we will use the properties of the 7 and a matrices 
and their identities. We will work in the representation of the 7 used by Bjorken and Drell [2J: 



" 1 " 


y = 


" a'' 




' 


1 " 


-1 


-a' 


7^ = 


1 






(41) 



and use the typical representation of the Pauli matrices: 



a 



' 


1 " 




" -r 




" 1 " 









i 




-1 


1 



(42) 



The Pauli matrices obey the identities 



0-^0-^ = 2 (j^ a^a^ = i cr^ = i 

(j^ (T^ = —i etc. 



(43) 



so therefore the 7 obey the identities: 



7^ 7^ = i 



-I a 



etc. 



(44) 



Here, the cr are 4 x 4 block-diagonal matrices formed from the 2x2 Pauli matrices. The 7 and a also 
have the properties: 

{a'f = l {jy = l {l'f = -l (45) 



Now consider the following: 



+ (7^)2 x'^d2 + 7^ 7^ x'^di + ... 



which by virtue of the identities (l44l) and (145|) becomes 



(7-^ x-')(7* 9j) = —x'^di + ia^ {x^d2 — x'^d^) + 
= — r c?/9r + a^Li/h 



so that we have established the identity 



d ■ 

(7-' x^){—ihY di) = ih r— i a^Li 

or 



which we will find convinient to recast as 



7° 7-' x^ 



I 



(46) 



When we write the 'j^ x^ object explicitly as 














1-2 

X X — I X 

x^ + ix"^ —x^ 



—x^ —x'^ + ix"^ 



1 • 2 

—X —IX 



X 








(47) 
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it becomes apparent that 

so that we can make the following statement about the inverse matrix 



7-' 



(48) 



Then, using (145!) with (1451) we obtain 



7° 7-^' x-' 



d .a' Li 

-in — + 1 

or r 



(49) 



so that the Dirac equation 
becomes 



(ihdo — mc 7°) ip = —ih'j^ Y di "ip 



This establishes Eq. ([1]), with the object e defined as 



-iTi — + 1 

or r 



7° 7'' x^ 



(50) 



When the explicit form of 7-^ x-' from ( 147[ 1 is combined with the definition of the spherical coordinates 

x^ = r sin 9 cos 
x"^ = r sin 9 sin 



X 



r cos 9 



We find that the explicit form of e becomes 







cos 9 sin9e~''' 

sin9e^'^ —cos 9 



cos 9 sin9e *" 
sin9e^'^ —cos 9 













(51) 



Note how in the definition of the e object, Eq. [50] and in other relations, the summation indicies j are 
both raised in the 'j^x^ object. In a simple flat space metric, this does not create a problem, so long as we 
note that since r = 'j-'xj, then e = —70^. 



Appendix B 

Here we will establish the properties of the 9^'^ 2-sp/nors when acted upon by A and o"*Lj. We will make 
use of the the properties of the angular momentum operators and spherical harmonics: 

L+Yi^rn = a/(/ - m)(/ + m + 1) 1^,^+1 /g2) 
L-Yi^m = ^/{l + m){l-m + l)Yi^m-i 
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La. — Z/, 



y /j + mj K,-i/2,r»j-l/2 

\/ j — rrij Yj_i/2,mj+i/2 



\/i + rnj(nij - 1/2) "K,- _ 1/2,™^. _ 1/2 + " '"j - '"jOO - 1/2 - "ij - 1/2 + 1) ^j- 1/2, -1/2 



Vi + m^ ^(j - 1/2 -m^ + 1/2)0 - 1/2 + m,- - 1/2 + 1) 1/2.™,. +1/2 " ("i^ +1/2)^3" 
(("ij - 1/2) + mj + - mj^a + mj)(j - m,j)) V-j_i/2,mj - 1/2 
(v^TTniJ y^Cj - m, )(j + mj) - (mj + 1/2) - ^j- 1/2, +1/2 

\/T+m~{mj - 1/2+ j - mj)Y,-_i/2,mj-i/2 



"3 '3-1/2, Tn,+l/2 



(53) 



y/j -rrijij + m^- - mj - l/2)Yj_^/2,mi+i/2 
h (j - 1/2) 91 



= h 



L+ -Lz J [ -y/j + mj + 1 yj+l/2,m-, +l/2 

y/j -rrtj +l( m,- - 1/2) y,-|-i/2,mj -1/2 - VT+ + 1 ^(7+ 1/2 + m,- + l/2)(j + 1/2 - mj - 1/2 + 1) y,-)-i/2,m,. -1/2 
^j-mj + 1 ^(j + 1/2 - +l/2)(j + l/2 + mj - 1/2 + 1) + i/2,r,ij +1/2 + ('"i + 1/2) ^j+ruj + 1 yj + i/2,mj-+l/2 
(("ij - 1/2) - + 1 - + rrij + 1^ (j + nij + l){i - + l)j ">Q + l/2,mj - 1/2 
(yj - + ly/(j - m,- + 1)0 + mj + 1) + (mj + 1/2)^^" + mj + l) V-^-+i/2,mj- + l/2 



y /j - rtij + I juij - 1/2 - j - m^- - l)yj+i/2,m,-i/2 
y/ITm~TT{j -mj + l+ rrij + l/2)Fj+i/2,m,+i/2 
-;i(j + 3/2)^7 



Now use the following identities to calculate the action of A: 



(54) 



cos e Yi^r, 



(2/+l)(2Z+3) 



(2;+i)(2;+3) 



{l-m+l){l+m+l) / (l-'m){l+m) 

-f/+l,m Y (2i-l)(2Z+l) - ■ 

" y 
(2i-l){2/+l) ^l-l,m+l 



(;-m)(/-m-l) 



(2i+l)(2Z+3) -f'+l,™-! 



(/+m){i+m-l) 
(2i-l)(2i+l) 



m— 1 



(55) 



cosQ sinO e 
sinO e***^ —cos9 



-^/j + rrij + l Yj+i/2,mj+i/2 



/ (j + 1/2 - mj-+ 1/2 + 1)0 + 1/2 + mj -1/2+1) /y + 1/2 -m,- + 1/2)0 + 1/2 + 1 
■111/ ^ ^3 + 3/2,™^- -1/2 + \ 



(2j+2)(2i+4) 



1/2) 



(2j)(2j+2) 



'j-1/2. 



",-1/2) 



+ 1/2 - mj - 1/2 + 1)0 + 1/2 - mj - 1/2 + 2) 



- + 1 

+ y'j + "ij + 1 

-mj- +1)^ (2j+2)f: 



(2i + 2)(2j+4) 



■^3+3/2, m,- -1/2 



(j + 1/2 + mj + 1/2)0 + 1/2 + mj + 1/2 - 1) 



(2j)(2j+2) 



^^3-1/2,™ .-1/2 



+ 1/2 + mj - 1/2 + 1)0 + 1/2 + mj - 1/2 + 2) 



(23 +2) (23 +4) 



i'3+3/2,mj+l/2 ■ 



+ 1/2 - mj + 1/2)0 + 1/2 - + 1/2 - 1) 



(23) (23 +2) 



^^3-l/2,m,+l/2 



0' + 1/2 - 



- 1/2 + 1)0 + 1/2 + mj + 1/2 + 1) 



1/2 

^jy^ + + "13 + 1) 



(I3 + 2)(2j + 4) 
/3 + 1/2+m,- 
V (23 + 2)6 



''3 + 3/2, ,„,j. +1/2 - 



+ 1/2 - mj - 1/2)0 + 1/2 + mj + 1/2) . 
(23)(23+2) 



'^3-l/2,m,-+l/2 I 



-1/2 1 



- ^ 



■j + 1)/ 



3 + l/2-m,--l/2 



(2;)(2;i2)"' + 0' + + 1)/ •'"(23'+2)?23) 



3 + l/2-m,--l/2 



1^3-1/2, m,- -1/2 



■^3-1/2, m,+l/2 



2i 



^ ^/T+m] Yj^i/2,mj-l/2 
^2^ - m^- K,_l/2,mj+l/2 



Because = 1 this then implies that 

H = \l-^(^7 (57) 
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